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Abstract 

The total perturbative contribution to the free-energy of hot SU (3) 
gauge theory is argued to lie significantly higher than the full result 
obtained by lattice simulations. This then suggests the existence of 
large non-perturbative corrections even at temperatures a few times 
above the critical temperature. Some speculations are then made on 
the nature and origin of the non-perturbative corrections. The anal- 
ysis is then carried out for quantum chromo dynamics, SU (N c ) gauge 
theories, and quantum electrodynamics, leading to a conjecture and 
one more speculation. 



1 Introduction 

The most convincing evidence for a phase transition in thermal Yang-Mills 
theories is provided by direct lattice simulation of the partition function, 

Z = Tre- pH , (1) 

where (3 = l/T is the inverse temperature. Over the years the lattice 
data for SU(3) theory, the purely gluonic sector of quantum chromodynam- 
ics (QCD), has become incresingly accurate, with various systematic errors 
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brought under control [[I]]. Figure (1) shows the normalized free-energy den- 
sity, F = —T\nZ/V, of SU(3) gauge theory taken from the first of Ref. 0. 
Plots such as this have supported a picture of a low-temperature phase of 
glueballs melting above some critical temperature to produce a deconfined 
phase of weakly interacting gluons: As the gluons are liberated the num- 
ber of degrees of freedom increases causing the free-energy density to rise, 
while asymptotic freedom guarantees the gluons are weakly interacting at 
sufficiently high temperature. 

Though the numerical data for QCD is less accurate, due to technical 
difficulties in simulating fermions, the accumulated data continues to support 
a phase transition. It is generally believed that this is a transition from 
a low-temperature hadronic phase to a high-temperature phase of quarks 
and gluons. This "quark-gluon plasma" is the new phase of matter which 
experiments at Brookhaven and CERN hope to detect in the near future. 

For the most part of this paper the focus will be on pure £77(3) theory, 
since the accurate lattice data allows a direct comparison with theory. Re- 
ferring to Fig.(l), there is one feature which is ignored by some, commented 
on by many and which has bothered a few. While there is little doubt that 
at infinite temperature a description in terms of gluons is tenable, this is less 
clear at temperatures a few times the critical temperature, T c ~ 270MeV. 
For example, at 3T C , the curve lies 20% below that of an ideal gas of gluons. 

What is the origin of this large deviation ? Is it due to (i) perturba- 
tive corrections to the ideal gas value, (ii) non-perturbative effects in the 
plasma, (iii) an equally important combination of (i) and (ii) or, (iv) is this 
an irrelevant question arising from an improper insistence of describing the 
high-temperature phase in terms of weakly coupled quasi-particles ? 

Several viewpoints have been expressed in the literature. Some believe 
that the deviation is mainly a non-perturbative correction to a gas of weakly 
coupled gluons and parametrize it in terms of a phenomenological "bag con- 
stant" . Others have attempted a phenomenological description of the high- 
temperature phase in terms of generalized quasi-particles. For a discussion 
and detailed references to these phenomenological approaches see, for exam- 
ple, 0. On the other hand, a few have suggested that the best consistent 
description of the high-temperature phase might be in terms of novel struc- 
turesSl. 

In order to help discern among the various possibilities, this paper will 
focus on estimating the total perturbative contribution to the free-energy. 
It is important to first agree on some terminology so as to avoid confusion 
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due to an overuse of some phrases in the literature. The partition function 
depends on the Yang-Mills coupling g, and has a natural representation as 
an Euclidean path-integral || 

Z(g) = J D<f> e -fo dT J d3x c ^ x ^ (2) 

where collectively denotes the gauge and ghost fields, and C the gauge-fixed 
lagrangian density of SU(3) gauge theory. An expansion of this path- integral, 
and hence the free-energy density, around g = leads to the usual Feynman 
perturbation theory and contributions of the form g n , with infrared effects 
occasionally generating logarithms multiplying the power terms, g n (ln(g)) m . 
These terms will be called perturbative. What is invisible in a diagramatic 
expansion about g = are terms like e~ x l g , associated for example with 
instantons M. Such terms, which are exponentially suppressed as g — ► 0, 
will be called non-perturbative. 

Note that at non-zero temperature, odd powers of the coupling, such as 
g 3 , appear ||. These are perfectly natural and represent collective effects in 
the plasma. Though they sample interesting long-distance physics, mathe- 
matically they fall into our definition of perturbative corrections. Similarly, 
Linde |J had shown that at order g 6 the free energy receives contributions 
from an infinite number of topologically distinct Feynman diagrams. Though 
the calculation of that contribution is difficult, it is possible in principle 0, 
and anyway does not qualify as a non-perturbative contribution according to 
the definition above. 

Following the heroic work of Arnold and Zhai, a completely analytical 
calculation of the free-energy of thermal gauge theory to order g 5 has been 
obtained f| 0|. For SU(3) gauge theory the result can be summarized as 
follows, 

£- -fe) + 30er + (6, 51 ng) +2 3, 2 - 2 0,3 ln (£)) 

where F = ~ 8 ^ r4 is the contribution of non-interacting gluons, a = g 2 /4n, 
and Jx is the renormalization scale in the MS scheme. Unfortunately @) is 
an oscillatory, non-convergent, series even for a as small as 0.2, which is close 
to the value of physical interest. A plot of (§) at different orders, at the scale 



p, = 2ttT is shown in Fig. (2). The poor convergence of (|3]) does not allow 
a direct comparison of the perturbative results with lattice data. Further- 
more, the result (0) is actually strongly dependent on the arbitrary value of 
ft. Inspired by the relative success of Pade' resummation in other areas of 
physics, Hatsuda |J and Kastening |I| studied the Pade' improvement of 
the divergent series (|^). Their conclusion was that the convergence could be 
improved, and the dependence on the scale fi reduced. However they did not 
attempt a direct comparison of their improved results with the lattice data, 
though Hatsuda did conclude that for the case of four fermion flavours, the 
deviation of the fifth-order Pade improved perturbative results from the ideal 
gas value was less than 10%. 

Not all seem to agree that a resummation of perturbative results as in 
H [HJ sheds sufficient light on the lattice data. For example, Andersen, et.al. 
Tl"| and Blaizot, et.al. [F2| have abandoned the expansion of the free-energy 



in terms of any formal parameter, but use instead gauge-invariance as the 
main guiding principle to sum select classes of diagrams. Their low order 
results seem to be close to the lattice data, but unfortunately because of the 
complexity of the calculations and the absence of an expansion parameter, 
it is not at all obvious what the magnitude of 'higher order' corrections is. 



A completely different approach has been taken by Kajantie, et.al. In [13 
an attempt has been made to numerically estimate the net contribution of 
long-distance effects, summarized in a dimensionally reduced effective theory, 
to the free-energy density. As will be discussed later, the calculation of [TB| 
probably contains some of the non-perturbative effects defined above but 
might miss out on some others. 

This author believes that the declared demise of information content in 
perturbative results such as @ is premature. In Ref.fl4|| a resummation 
scheme was introduced to obtain an estimate of the total perturbative con- 
tribution to the free-energy density of SU (3) theory. The methodology of 



Ref. |b| has been further developed and applied to other problems in |L5], [L6 



is 



In Sees. (2-5), an explicit and improved discussion of the results in |14 
given, leading to the conclusion that the total perturbative contribution to 
the free-energy density lies significantly above the full lattice data. In Sec. (6) 
I discuss the consequent magnitude of non-perturbative contributions, and 
speculate on their possible origin. Sec. (7) contains an analysis of the per- 
turbative free-energy of generalized QCD, with Nf fermions, and a brief 
comparison with the available lattice data which is less precise. SU(N C ) 
gauge theory is discussed in Sec. (8), and an apparently universal relation 
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noted. Sec. (9) considers quantum electrodynamics (QED) and some specu- 
lations about its high-temperature phase. A summary and the conclusion is 
in Sec. (9). 



2 The Resummation Scheme 

The truncated perturbative expansion of the normalized free-energy density 
can be written as 

TV 

S N (X) = l + J2fnX n , (4) 

n=l 

where A = (01/77)2 is the coupling constant, and where following [10|, 
possible logarithms of the coupling constant are absorbed into the coefficients 
f n . The poor convergence of (f|) is obviously due to the large coefficients at 
high order. Indeed such divergence of perturbative expansions is generic in 
quantum field theory and one expects the coefficients /„ to grow as n\ for 



large n [[L7]]. This leads to the introduction of the Borel transform 



N f n 

B N {z) = 1 + E ^ (5) 



n=l U - 



which has better convergence properties than (f|). The series (f|) may then 
be recovered using the Borel integral 

1 r°° 

S N (X) = - dz e~ z ' x B N (z) . (6) 



A Jo 

The logic of Borel resummation is to define the total sum 5(A) of the per- 
turbation expansion as the N — > 00 limit of @. This of course requires 
knowledge of B (z) = B^z) and the existence of the Borel integral. Lacking 
knowledge of the exact B(z), one therefore attempts to reconstruct an ap- 
proximation to 5 by replacing the partial series B^(z) in @ by a possible 
analytical continuation thereof. A simple and popular method to achieve this 
is to use Pade' approximants, leading therefore to a Borel-Pade' resumma- 
tion of the series (|4]). This method will be briefly discussed in Sec. (5). Here 
instead I will proceed as suggested by Loeffel |L8| and change variables in (||) 
through a conformal map. For a positive parameter p, define 

y/l+pZ-l , , 

w \ z ) — n 7 {<) 
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which maps the complex z-plane (Borel plane) to a unit circle in the w-plane. 
The inverse of (0) is given by 



Aw 
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The idea |19| is to rewrite 



p (1 — w) 2 
in terms of the variable w. 



S N (X) 



z n is expanded to order N in w and substituted into 

(2n + k - 1] 



(8) 

Therefore, using 
The result is 




k\ (2n - 1) 



e- z/x w(z) ik+n) dz , (9) 



where w(z) is given by (|7|). Equation (Qj represents a highly nontrivial re- 
summation of the original series @ |15|] . In the pioneering application of the 
Borel-conformal-map technique in condensed matter physics by Le-Guillou 
and Zinn- Justin |19| , the paramter p was a fixed constant which determined 
the precise location of the instanton singularity at z = — 1/p. In some more 



recent QCD applications p0| , the fixed constant p determines the ultraviolet 
renormalon singularity closest to the origin in the Borel plane PT |. 



The novelty introduced in Ref.fl4fl and futher developed in |15] was to 
consider p as a variational parameter determined according to the condition 

dS N (X ,p) 



dp 







(10) 



where Ao is a convenient reference value. For the problem at hand, because 
/2 is negative, the solutions p(N) to ( |I0|) will be postions of minima |15| . 
Denote the value of (||) at p = p(N) by Sat (A). Notice that although p(N) is 
determined at the reference value Ao, S'jv(A) is defined for all A. The reason 
why this is sufficient has been explained in [|l5j and will be discussed further 
in the next section. 

It must be stressed that, in general (see |15|]), the variational parameter 
p(N) is not related to possible singularities of B(z). Rather, it is determined 
according to the extremum condition (0). Thus the presentation here is a 
slight departure from that in [T3j and represents the developments in Ref . [T3 



Sufficient information is now available to construct the resummed approx- 
imants SV from iV = 3 up to iV = 5 in the next section. The approximant Si 
of course does not exist since fi = 0, while S2 cannot be constructed because 
no solution exists to equations such as (|10|) at the first nontrivial order. For 
N > 2, the solutions p(N) will be positions of global (local) minima if the 
sign of /tv is positive (negative) [ O . 
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3 SU(3) : Resummation up to Fifth Order 



In order to make contact with lattice data which show a temperature depen- 
dent curve, one must use in (f|) a temperature dependent coupling. Let us 
begin by using the one- loop running coupling defined by [ITTl p~2]j, 



A(c,x) = ^^_, (11) 
HLfc, x) 



where L(c, x) = ln((2.287rcx) 2 ), c = fi/2nT and x = T/T c , with T c ~ 
270MeV the critical temperature which separates the low and high tem- 
perature phases [1], [2^] . 

Fixing first the reference values Co = 1, Xq = 1, which fixes the reference 
value of A , the results of ([T(J are: p(3) = 3.2, p(4) = 7.9, p(5) = 13.7. The 
curves for «Sat(A) are shown in Fig. (3a) at the renormalization scale c = 1. 
Notice the behaviour S5 > S 4 > S3 and how these all lie significantly above 
the lattice curve in Fig.(l). The results do depend on the renormalization 
scale, denoted here by the dimensionless parameter c. It has been suggested 
fTT| , P~2j1 that a suitable choice for such a parameter is 0.5 < c < 2, corre- 
sponding to 7rT < /2 < AttT. Certainly this is the natural energy range for 
the high-temperature phase. Figure (3b) shows the mild dependence of S5 
on the renormalization scale. 

The results above were obtained by solving (|It]) at the point Co = xq = 1. 
Now consider changing the reference values to c = 1, x = 3, that is, a more 
central value for the temperature. The solutions are: p(3) = 3.2, p(4) = 
7.8, p(5) = 13.4. These values are hardly different from those above. This 
is firstly due to the fact that (§) is a much slower varying function of the 
coupling than the original divergent series. Furthermore, for the present 
problem, the coupling itself varies slower than logarithmically with c and x 
(the c and x dependence of the coefficients / 4 and f 5 is also only logarithmic). 
The curves for the re-optimized Sn are essentially identical to those shown in 
Fig.(3a,b), the difference being only at the fifth decimal point. For example, 
the value of S$ in Fig. (3a) at x = 3 is 0.938684, while that for S5 optimized 
at Xq = 3 (and hence evaluated at p(5) = 13.4), is 0.938672 at x = 3. This 
confirms the assertion in [fbl that the results are quite insensitive to the 



exact reference values chosen to solve (|I0|). 

We now proceed to test the sensitivity of the results to the approxima- 
tion used for the running coupling (111]). The approximate two-loop running 
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coupling is given by |TI, VI 



Mc x) - 2 (i - 51 , 12) 

( ,_ 121 l(c,x) ) (12) 



with the symbols having the same meaning as before. In Fig. (4), the one- 
loop running coupling (|TT| ) and the approximate two-loop running coupling 
([12]) are plotted at c = 1 to show their difference. At a; = 3, the value 
for the approximate two-loop coupling is about 20% lower than the one- 
loop result. Nevertheless, because of the above-mentioned property of the 
resummed series, we shall see that the final results to do not shift dramati- 
cally. Using (0), the solution of (|i~0| ) at the reference point Co = l,£o = 3 
are: p(3) = 3.2, p(A) = 7.6, p(5) = 13.1. The corresponding curves shown 
in Fig. (5a) have moved up slightly compared to those in Fig. (3a). The "two- 
loop" value of S 5 (c = l,x = 3) = 0.9473 should be compared to the "one- 
loop" value 0.9387 obtained above. The mild renormalization scale depen- 
dence of the new 5*5 is shown in Fig. (5b). 

In summary, it has been demonstrated that the resummed approximants 
S3, S4, S*5 all lie significantly above the lattice data and satisfy the mono- 
tonicity condition S5 > S4 > S3. The result is insensitive to the reference 
value used to solve (|I0D, for the range of interest 0.5 <c<2, 1 < a; < 5. The 
result is also insensitive to the approximation used for the running coupling 
constant and in fact better approximations for the coupling seem to move 
the values of Sn further away from the lattice data. Finally it should be 
noted that the values Sn also appear to converge as N increases. 

The only way to force the values of S^ down closer to the lattice data 
is to choose very low values for the renormalization scale, c ~ 0.05. Of 
course this is not only unnatural but increases the effective value of the 
coupling contant beyond what one would believe is physically reasonable for 
a perturbative treatment. That is by making an artificially low choice for 
the renormalization scale, one cannot escape the conclusion stated in the 
abstract of large non-perturbative corrections ! 



4 Higher Order Corrections 

Due to technical complications the sixth order contribution, A 6 , to the free- 
energy density has not been calculated although an algorithm for it exists 
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0. There is a misconception that because that contribution is due to an 
infinite number of topologically distinct diagrams, its value must be very 
large. A counter-example is provided by the magnetic screening mass ||, 
which suffers from the same disease but whose approximate calculations in 
the literature show it to be of ordinary magnitude f23 |. 

Having said that, let us see what is the worst that can happen. It has 
been suggested [|24j that Pade' approximants can be used to estimate the 
next term of a truncated perturbation series. That is, after approximating 
the truncated series by the ratio of two polynomials, the Pade' approximant 
is re-expanded as a power series to estimate the next term in the series. 
Well, why not also use Borel-Pade' approximants for the same purpose ? 
Using the fifth order result (|3|) together with the two-loop running coupling 
(12|), and choosing the central values c = 1, x = 3, all fifth order Pade and 
Borel-Pade approximants were constructed and then re-expanded to give an 
estimate of the coefficient fa. The largest value obtained was 30, 000 and the 
smallest —30, 000. Note that the fifth order coefficient at c = 1 is —800, so 
the estimated magnitude of f 6 is about 37 times larger. Since the coupling 
A is about 0.2 at x — 3, c = 1, the total value of the sixth order contribution 
to (|3|) is therefore estimated to be almost 8 times in magnitude compared to 
the fifth order contribution. These are big numbers and should be expected 
to cause some damage. 

Using (Q) with / 6 = 30,000, the two-loop coupling (0), and solving 
flTT)| ) at the reference point c = 1, x = 3 gives p(6) = 19.75 and S e (c = 

1, x = 3) = 0.9490. Repeating for / 6 = -30,000 gives p(6) = 19.5 and 
Sq(c — 1, x — 3) = 0.9489. Notice the negligible change in the value of 
5*6 even when wildly differing values have been used for / 6 . Those values 
should be compared with the fifth order approximant of Fig. (3a), which gives 
5*5(0 = l,x — 3) = 0.9473. The large estimated sixth order corrections to 
the divergent perturbation expansion (||) cause a change of only 0.002 to the 
values of the resummed series, and more importantly the shift is upwards, 
S$ > S5, preserving the lower order monotonicity, regardless of the sign of 



Kajantie, et.al. [13| have suggested that the sixth order contribution, 
f$\ 6 be of order 10. For a coupling A ~ 0.2, this translates into the astro- 
nomical value ±156250 for / 6 . Solving (|10|) at Co = l,a;o=3 gives p(6) = 19.8 
and Sq(c — l,x — 3) = 0.9491 for the positive fa, and p(6) = 19.4 and 
S 6 (c — l,x — 3) = 0.9489 for the negative f 6 . Despite the anomalously large 
value of the sixth order contribution proposed in |13[], the conclusion here is 
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still S 6 > S 5 , and an increment of only 0.002. 

To highlight the above result in a more dramatic way, suppose the sixth 
order coefficient vanishes, fe — 0. Then because of the non-trivial way the 
resummation is done in (^J, the solution to (|10D for N = Q will still be 
different from that of iV = 5. At Cq — 1, Xq — 3 I find p(6) = 19.5 and then 
S 6 (c — 1, x — 3) = 0.9489 at four decimal places, which is almost identical 
to the values obtained above for various large values of /g. This sounds 
incredible but is actually not once one remembers that large corrections to the 
divergent series @ do not translate into large corrections to the resummed 
series (||). In fact those large values are suppressed in various ways. Firstly, 
in the re-organization of the series in flj^), less weight is given to higher order 
corrections. Secondly, the variational procedure chooses values of p(N) which 
in this example increase with N, and so suppress further the value of Sn- 

More understanding of the above results can be obtained through a large 
iV analysis carried out for the general Eqs. (|9],|T0|) in |15|]. It was shown in 



T5| that if p(N) increases for the first few values of N, then that trend will 
continue. Let c(N) = p(N + l)/p(N). In the large N and large p(N) limit 
one can show that [15| 

1 11 

1 — Z7at\ + Znvrv? ■ (13) 



c(N + l) c(N) c(N)< 

One consequence of this is that c(N + 1) < c(N) and c(N) — > 1 + as N — > oo. 
This is indeed observed for the present problem already at low N. Numeri- 
cally, ([131) too is not a bad approximation at small N. In fact using the val- 
ues found for p(N) in the last section, one has p(4)/p(3) = 7.6/3.2 = 2.375. 
With this as input for c(3), (|1^) gives the estimate c(4) ~ 1.32234, to be 
compared with the actual value p(5)/p(4) = 13.1/7.6 = 1.7. Next using 
c(4) = 1.7 as the exact input, (pf) gives c(5) ~ 1.32201 and thus an estimate 
of p(6) ~ 1.3 x 13.1 = 17.03. On the other hand, using p(6) ~ 19.5, as de- 
termined by various estimates above, gives c(5) = 19.5/13.1 = 1.49 and then 
through d|) the estimate c(6) ~ 1.28 and hence p(7) ~ 1.28 x 19.5 = 25. 
The Eq.(|13|) was derived in [|15] for the case fi ^ 0. For the present case 



where fi = one will actually obtain the slightly more accurate equation 

1 1 1 

l-77^ + -7^> (14) 



c(iV + 1) 2 c(N) c(N) 3 

but in the large N limit where c — > 1 + this is clearly equivalent to ( jT3[) . 
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Note that the recursion relations (|TB| , TA) make no explicit reference to 



the values of the f n which in the derivation in [|1J were assumed to be 
generic, that is, diverging at most factorially with n. Indeed that fact that 
various different assumptions about the value of f$ earlier in this section led 
to essentially the same value for p(6) ~ 19.5 supports the f n independence 
of (PD already at iV ~ 6. 



From the general analysis in [fig] , one also deduces that for large N and 
large p(N), the monotonicity ASn = SV+i — Sn > is guaranteed by the 
fact /2 < 0, and that ASn ~ 1/iV 3 as N — > oo. Since the explicit N < 5 
calculations and the estimated N = 6 result already support c(N) > 1 and 
large values of p(N) at low N, this suggests that the continued monotonicity 
and rapid convergence of the Sn is assured by the large N analysis. 



5 Lower Bound and Other Estimates 

From the explicit low iV calculations, and the large N analysis, one concludes 
that for N > 2, 

Sn < Sjv+i (15) 

for all N, and furthermore, the difference Sn+i~ Sn decreases as N increases, 
showing a rapid convergnce of the approximants. However, in general, it is 
not quite correct to say that the approximants converge to the total sum of 
the series [0. For each N, let p*(N) be the value of p that is optimal, that 
is, it is the value which when used in (Dp gives the best estimate of S, the 
total sum of the series. Define, S N = Sn(\p*(N)). Then for those p(N) 
which are positions of global minima one has by definition, 

S N < S N (16) 

It is S N which presumably converges to S as N — > oo. (This implicitly 
assumes that the sub-sequence of global minima is infinite: That is, given 
any positive integer iVo, there is some n > N for which /„ is positive.) 

Hence if one accepts the two assumptions above, then combining flT5D 
with flU), 

S N < S (17) 

for all iV > 2, and one may conclude that the Sn are lower bounds to the 
sum of the full perturbation series. 
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In particular that conclusion implies that the N = 5 curve in Fig. (5a) is a 
lower bound on the total perturbative contribution to the free-energy density 
of hot SU(3) theory. The statement has three qualifications. Firstly, it 
involves the technical assumptions mentioned above. Secondly, as discussed 
before, better approximations to the running coupling can move the bounds, 
but it was seen that a 20% improvement in the coupling shifted the bound 
upwards by 1%. Thirdly, the bounds shift by ±1% when the renormalization 
scale is varied by a factor of two from its central value p, = 2-kT. Thus 
it might be more appropriate to call the bounds as "plausible soft lower 
bounds" with an uncertainty depicted in Fig. (5b). 

Given that the lower bound obtained above involves some unproved tech- 
nical assumptions, it is useful to compare the above results with those ob- 
tained using different resummation schemes for the divergent series (|]). I 
briefly state here the main results obtained using a Borel-Pade' resumma- 
tion of (U), with the two- loop approximation for the coupling (|T2"D and the 
central value c = 1. The approximants will be denoted as [P, Q], referring 
of course to the particular Pade' approximant used for the partial Borel se- 
ries (|D constructed from (|J). The only approximants which did not develop 
poles and which gave a resummed value below one in the temperature range 
2 < x < 5 were [1,2], [2,1] and [2,2]. These are displayed in Fig. (6). The 
[3, 2] and [4, 1] approximants did not develop poles but gave a value above 
one. If the approximants which developed a pole are defined through a princi- 
pal value prescription, then the lowest value was given by [2, 3]: 0.91 — > 0.94 
as x increased from 2^5. The [1,3] and [1,4] approximants gave values 
above 0.98 in the range of interest while [3, 1] gave a value above one. 

Thus in the Borel-Pade' method, the minimum estimate for the fifth order 
resummed series is given by the principal value regulated [2,3]. The highest 
values were all above one. If one keeps only the fifth order estimates below 
one (thus giving a very conservative lower value), then the average of the 
[2,3] and [1,4] is greater than 0.94 for the entire range 2 < x < 5. At 
x = 3 the estimates are 0.95 ± 0.04. Of course including also the values 
above one would push this average higher. Clearly the Borel-Pade' estimates 
are comparable to the bounds obtained using the resummation technique of 
Sec. (2) and should reassure some readers about the novel resummation used 
here. 

For completeness, I mention an alternative way of thinking about diver- 
gent series such as (f|). For QCD at zero temperature, a paradox is that one- 
loop results give remarkable agreement with experimental data even when the 
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energy scale is relatively low. As the running coupling is then large it is not 
obvious why higher-loop perturbative corrections are suppressed. It has been 
suggested |25], |26| that the explanation might lie in the probabale asymptotic 
nature of the QCD perturbation series. Recall that in an asymptotic series 
the best estimate of the full sum, at a given value of the coupling, is obtained 
when only an optimal number of terms is kept and the rest discarded (even 
if they are large). Thus if one knew the general behaviour, at least at large 
order, of the series and assumed that it was asymptotic, then one could 
have obtained a reasonable estimate of the full sum by simply adding the 
optimal first few terms. What has been done in the previous sections, and 
this is what various resummation schemes try to do, is to instead sum up the 
whole series to get an even better estimate of the total perturbation series 
(and this has the greater advantage of giving a good result for a large range 
of couplings). Also note that thinking of (f|) as an asymptotic series does not 



say anything about explicit non-perturbative corrections [25, 26 



6 Non-Perturbative Corrections 

The total perturbative contribution to the free-energy density of SU(3) gauge 
theory has been argued to be close to, or above, the N = 5 curve in Fig. (5a). 
A residual uncertainty that could lower the curve of Fig. (5a) is the exact 
value of the renormalization scale. For a natural range of parameters, the 
lower curve in Fig. (5b) is the result. On the other hand the full result as given 
by lattice simulations is shown in Fig.(l). Lattice errors have been stated to 
be under 5% [[![]. Taken together, the conclusion appears inescapable : Even 
at temperatures a few times above the transition temperature, there are large 
negative non-perturbative contributions to the free- energy density. For exam- 
ple, at T = 3T C ~ 700MeV, the lattice results for the normalised free-energy 
density are 0.8 ± 0.04 while the lower bound on the perturbative contribu- 
tion is 0.947 ± 0.007, implying a minimum non-perturbative correction of 
10% (and as high as 20%). 

Thus an answer has been given to the questions raised in the introduc- 
tion. The deviation of the lattice data from the ideal gas value is apparently 
caused mainly by non-perturbative corrections, with perturbative corrections 
accounting for a much smaller amount. At T ~ 3T C the relative contributions 
are ~ 15% and ~ 5%. 

I speculate now on possible sources of the non-perturbative corrections. 
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Firstly there are the familiar instantons, already present in the classical ac- 
tion, and which contribute terms of the order e _1//A . Secondly there are the 
magnetic monopoles. There is by now overwhelming evidence that confine- 
ment at zero temperature is caused by the t'Hooft-Mandelstam mechanism 
of condensing monopoles (the dual superconducting vacuum). Thus it is 
possible that the monopole condensate has not completely melted above the 
critical temperature. Note that since the classical theory does not support 
finite energy monopoles, these must be of quantum origin, and so their con- 
tribution might be larger than those of the instantons. 

In fact contributions which are exponentially small but much larger than 
those of the instantons are suggested by the Borel resummation itself. It is 
known that Yang- Mills theories are not Borel summable |2T]|. That is, 



the function B{z) has singularities for positive z, making the Borel integral 
ill-defined. One can nevertheless define the sum of the perturbation series 
using the Borel integral if a prescription is used to handle the singularities. 
It is generally believed that the prescription dependent ambiguity disappears 
when explicit non-perturbative contributions are taken into account for the 
physical quantity in question. Indeed the nature of singularity itself suggests 
the form for the non-perturbative contribution. If there is a pole at z = q, 
then the non-perturbative contribution will be of the form ~ e~ g ^ x , which is 
larger than the instanton contribution for small A. An explicit mathematical 
model which illustrates the interplay between Borel non-summability and 
non-perturbative contributions has been given in |15[| . 

Notice that the non-perturbative corrections suggested by the Borel method 
at non-zero temperature are very different from those at zero temperature. 
In the latter case the expansion parameter is g 2 and so the contribution is 
~ e~~ q / 9 , which translates into a power suppressed contribution ~ l/(Q) b 
when g 2 is replaced by the running coupling ~ l/ln(Q/A). In cases where 
the physical quantity can also be analysed using the operator-product ex- 
pansion (OPE), these power suppressed contributions to perturbative results 



correspond in the OPE picture to vacuum condensates |20|, . 

At nonzero temperature, since the natural expansion parameter is A = 
\J g 2 /A-k 2 ~ 1/ Jla(T/A), one does not get a simple power suppression from 

e -1 A. Nevertheless, the analogy with zero-temperature results suggests that 
such contributions might be due to some condensates. Thus the conven- 



tional condensates discussed for example in |L3| are plausibly part of the 



non-perturbative contributions. The form however suggests even more novel 
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condensates. These might be, for example, those of DeTar [[| or Pisarski 

It is worth noting an explicit instance of a theory displaying exponentially 
small non-perturbative effects which are larger than those due to standard 
solitons. In fundamental string theory where the coupling is g, there are 
the usual solitons but there are also novel " D-instantons" which give a larger 
contribution e -1 / 9 [27]. I also mention in passing the recurrent and intriguing 



relationship between gauge theories and strings p8| which leads one to won- 
der whether that is a possible route to understanding the non-perturbative 
structure of hot gauge theories. 

Using the N = 5 curve of Fig. (5a) as a reasonable estimate of the full 
perturbative result, and assuming a nonperturbative component of the form 

S np =je-«l\ (18) 
as suggested by the Borel method, one can determine the constants A and q 



by comparing the lattice data of Fig.(l) with the perturbative result. In [14 
it was shown that, 

Slatt = S pert - -L e8.7-2.62/AM (ig) 

A(x) 

where Si a u represents the lattice data for the free-energy, and S pert the re- 
summed perturbative result, both normalized with respect to the ideal gas 
value, and X(x) is given by ( |Ti2"D at c = 1. Eq.(|TjJ) is a phenomenological equa- 
tion of state for the free-energy which generalises the usual discussions in the 
literature where the second term on the right-hand-side of ([19|) is called a 'bag 
constant'. In this case the 'constant' is really temperature dependent and 
represents a non-perturbative contribution to the free-energy that vanishes 
at infinite temperature. It is important to note that the non-perturbative 
contribution is negative, since the perturbative result is above the full lattice 
data, and thus consistent with the usual interpretation in the literature. 

So far the discussion has implicitly assumed an additive picture of per- 
turbative and non-perturbative contributions, with both components clearly 
distinguished. It might be that in reality the best description of the high- 
temperature phase is in terms of completely novel structures |3], ||. In that 
case a forced expansion of those alternatives about A = must give some- 
thing like (|19D and the subsequent distinction between perturbative and non- 
perturbative contributions. The mathematical toy model of [|HJ illustrates 
this. 
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Within the framework of this paper, one can distinguish three versions of 
the popular concept of " quasi-particle" . Firstly there are the "perturbative 
quasi-particles" which are deformations of the gluon formed by a particular 
reorganisation of the perturbative Feynman diagram expansion. The results 
above suggest that if all the contributions of such quasi-particles to the free- 
energy are added up, the net result will lie above the lattice data, and only 
truly non-perturbative contributions, as defined in Sec.(l), may give the final 
agreement. Secondly there are the "nonperturbative quasi-particles", which 
are excitations about the nontrivial thermal vacuum that includes conden- 
sates, and so forth. Currently there is not sufficient control over the theory to 
construct these objects. Finally there are the " phenomenological quasipar- 
ticles" which simply aim to give numerical agreement with the lattice data 
within a simple ansatz. The ultimate justification for these phenomo logical 
constructs must surely come from the " nonperturbative quasiparticles" . 



7 QCD 

In this section lower bounds (within assumptions similar to those made pre- 
viously) are obtained for the perturbative free-energy density of hot SU(3) 
coupled to Nf flavours of fundamental fermions. As the essential features are 
very similar to the pure gauge case, the discussion here will be brief. The 
fifth order perturbative results in the MS scheme can be read off from the 
landmark papers ||. The approximate two- loop coupling that is used here 
is given by 

X{c,x) = = 1 - 2 — (20) 

^/A7i 2 p L(c, x) \ 2/3 2 L(c,x) J 

with 

llN c -2N f , s 

A) = — h^r^> ( 21 ) 



48tt 

A = ^y(^ N c-^ N cNf + SN f /N c ). (22) 



Following [0 I also assume a relative Nf independence of L(c, x), and thus 



use for it the same expression as used in ([11]). The extremization condition 



(|T0| ) is solved at the reference point cq = 1, xq = 3 for iV c = 3 and 1 < 
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Nf < 6. The results obtained are all extremely similar: In each case the 
convergence of the Sn is monotonic and rapid as in the pure glue case in 
Fig(5a). For this reason only the N = 5 curves are displayed in Fig. (7) for 
the various number of flavours. For comparison the pure glue result (Nf = 0) 
is also included. (Each curve has been normalized with respect to the ideal 
gas value for that number of flavours.) 

The N = 5 curves in Fig. (7) can be taken as plausible lower bounds, or 
estimates, to the total perturbative free energy density of QCD with Nf fun- 
damental fermions. Lattice results for QCD contain large systematic errors 
compared to those for SU (3) and so a precise comaprison is not possible. Af- 
ter making some assumptions about the size of the systematic errors, the au- 
thors in Ref. [2^| determine that for Nf = 2 the free-energy density lies about 
15 — 20% below the ideal gas limit. This is similar to the case of pure SU(3). 
Comparing this lattice estimate with the estimate on the perturbative result 
in Fig. (7) one is again led to suggest that there are large non-perturbative 
corrections to the naive picture of a weakly interacting quark-gluon plasma. 

Of course, given the physical relevance of QCD, it would be preferable to 
have more precise numbers from the lattice, and especially for other values 
of Nf. However it seems that a non-perturbative component of 10 — 15% at 
temperatures a few times T c is likely to be generic. 

When the draft of this paper was complete, I came across [22] which 
gives for Nf = 2 a T c /Ajj^ ~ 0.5, a factor of two lower than that for the 
pure gauge theory. This has the consequence that c in L(c,x) should be 
replaced by about c/2. However as the reader can surmise by now, this has 
hardly any impact on the results above, for this is equivalent to shifting the 
renormalization scale c by a factor of two, which as we have seen causes only 
a 1% shift of the curves. In any case this serves to remind that the lattice 
results for Nf are in a state of flux. 



8 SU(N) ( 

Define 

\<N) =(■ 

3 J \ti, 

Then the free-energy density of pure SU(N C ) theory up to fifth order |§. [7j 
is given by the expression (|3|), with (a/n) 1 ^ 2 replaced everywhere (including 



fN \ 1/2 /a\ 1 / 2 

W) = (ir) (-) • (23) 
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inside the logs) by X(N C ). Thus there is no explicit N c dependence of the 
free-energy density when written in terms of X(N C ). 

To examine the N c dependence of the new coupling fl23|), consider the 
approximate two-loop running coupling given by 



with 



/2r'irT\ 

L(T,A) = 21n^ , (25) 

& = W' (26) 

A = ^y^ N c)i (27) 



and with A = A(N C ) the SU(N C ) gauge theory scale parameter in the MS 
scheme. In L, the constant c' is Ji/2nT. Comparing the various equa- 
tions, one comes to the remarkable conclusion that the new running coupling 
X(N C ,T) will be independent of N c if the MS scale parameter A(N C ) is itself 
independent of N c when expressed in terms of some physical length scale. 
By comparing some data for iV c = 2, 3, 4, Teper ]30f has concluded that this 
is indeed the case. 



Therefore, accepting the result of [p0| , one deduces that the N = 5 curve 



in Fig. (5a) is a plausible lower bound, or estimate, to the total perturbative 
free energy density of hot SU(N C ) theory when the x-axis is interpreted as 
T/A instead of T/T c . This is then a universal relation (at least for low 
N c ), and one suspects that the corresponding full lattice results might also 
obey a universal curve, thus leading to the guess that the non-perturbative 
component of an SU (N c ) plasma is about 10 — 15% for temperatures a few 
times above the critical temperature. 



9 QED 

Though the fine structure constant a of QED is small at everyday energies, it 
is interesting to consider super-high temperatures where it will be large. The 
free-energy density of massless QED at temperature (T), has been computed 
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up to order a 5//2 in Refs. pH|, [| |7], Denoting as usual A = (a/rr) 1 ^ 2 , the 
normalised free-energy density at the MS renormalisation scale p, = 2tcT, is 
given by |JT], || 



F/F = 1 - 1.13636A 2 + 2.09946A 3 + 0.488875A 4 - 6.34112A 5 . (28) 

where F Q = ll7r 2 T 4 /180 is the free-energy density of a non-interacting plasma 
of electrons, positrons and photons. Figure (8a) shows the plot of ( P5| ) at 
different orders. The series diverges at large coupling (super-high tempera- 
tures), exhibiting a behaviour similar to that of Yang-Mills theory at low- 
temperatures. The convergence at large coupling can be improved by using 
the resummation technique (0-0). Using the coefficients from (p8|), the solu- 
tions of (|IUD at the reference value A = 0.5 are (minima): p(3) = 0.7, p(4) = 
1.75, p(5) = 3. 

The resummed series, with its much improved convergence, is shown in 
Fig. (8b). The N = 5 curve can be taken as a lower bound to the full perturba- 
tive result. If one assumes that the potential non-perturbative contributions 
lower the perturbative result, as happens in QCD, or are very small in mag- 
nitude, then one may conclude from Fig. (8b) that super- hot QED undergoes 
a phase transition. This speculated high-temperaure phase of QED might 
then be analogous to the low-temperature phase of QCD with various bound 
states. Or, it might resemble the alternative picture of low-energy QCD: that 
of flux-tubes [Q. It is unfortunate that no lattice or other non-perturbative 
information is currently available about the high-temperature phase of QED. 



10 Conclusion 

The phrase 'non-perturbative' is used often and loosely with regard to field 
theories at non-zero temperature. This has caused a great deal of semantic 
confusion and misunderstanding. For the purpose of uncovering the cause 
of the deviation of the result in Fig.(l) from the ideal gas value, it has been 
proposed to term 'perturbative' all power like (modulo logarithms) contribu- 
tions to the free-energy density. Such perturbative contributions follow from 
the usual Feynman diagram expansion of the parition function around zero 
coupling. 

For SU(3) gauge theory a plausible lower bound was obtained on the 
totality of such perturbative contributions to the free-energy density. The 
derivation of that lower bound using the variational conformal map involved 
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some technical assumptions, and so one may instead wish to consider it only 
as an estimate of the total perturbative contribution. The estimate is compa- 
rable to that obtained using Pade' or Borel-Pade' resummation methods and 
lies significantly higher than the full lattice result, thus suggesting that large 
and truly non-perturbative corrections exist. As discussed in Sec. (6), these 
non-perturbative corrections might include the usual instantons, magnetic 
monopoles, the usual condensates, and perhaps also more novel condensates 
and extended structures as suggested by the Borel form e _1 ^ A of the non- 
perturbative contributions. As to which of these possibilities dominates is 
an interesting question left for future work. 

The equation of state for hot SU(3) can be summarized by the phe- 
nomenological relation WA 



S latt = S pert - e 8 ' 7 - 2 - 62 ^) , (29) 

where Si att represents the lattice data for the free-energy, and S pert the re- 
summed perturbative result, both normalized with respect to the ideal gas 
value, and where A is the temperature dependent coupling (|l"2"D at c = 
1. There is a slight ambiguity in the estimate of the magnitude of non- 
perturbative corrections coming from the residual renormalization scale am- 
biguity of the resummed perturbative results. For the natural range irT < 
Ji < AttT, the ambiguity is less than one percent. Such an ambiguity be- 
tween the perturbative and non-perturbative components is understandable, 
as only the full physical quantity can be demanded to be scale indepen- 
dent, and not separately its perturbative and non-perturbative components, 
though the latter simplifying assumption is often made. 

By choosing an anomalously low value for the scale fx one can fit the 
purely perturbative results to the lattice data but at the price of a very 
large effective coupling constant, and thus an unsuppresed contribution from 
common place objects like instantons. Therefore this route does not offer an 
escape from the conclusion of large non-perturbative corrections. 

The results of this paper do not imply that the description of the free- 
energy density must be as in fl29|). Rather, completely novel descriptions 
in terms of various extended structures are possible |3|, M, but in a forced 
expansion of those alternative descriptions about g = 0, one must recover 



something like (|19|). On the other hand, the results here do imply that 
an accounting of the lattice free-energy density in terms of a subset of Feyn- 
man diagrams of the perturbative thermal vacuum might miss some essential 
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physics. 

Being restricted to an analysis of the free-energy density, the results 
here of course do not imply that every observable must have large non- 
perturbative contributions, but it is likely that this is true for most of the 
bulk thermodynamic quantities. 

An analysis was also carried out for generalized QCD with Nf quarks. 
Currently the lattice data is only approximate and so the conclusions are 
less definitive. However accepting the estimates in [p9[] , the conclusion is the 
same as before: There appear to be large non-perturbative corrections to the 
free-energy density of hot SU(3) gauge theory coupled to Nf quarks. Thus 
in particular, "quark-gluon plasma" seems to be an incomplete description 
even at temperatures several times above the transition temperature. 

A simple relation was noted for the perturbative free-energy density of 
SU(N C ) gauge theory in Sec. (8). That result, combined with the results 
for generalised QCD in Sec. (7), and the available lattice data, leads one to a 
universality conjecture: For the high-temperature (T > 2T C ) phase of SU(N C ) 
gauge theory coupled to Nf fundamental quarks, and for all moderate values 
of N c and Nf, at most 5% of the deviation of the free-energy density from the 
ideal gas value is due to perturbative effects while non-perturbative effects 
contribute a larger 10 — 15%. 

It has been speculated in Sec. (9) that QED might have an interesting 
high-temperature phase. 

Finally, the methods of this paper might be of some use for the study 
of supersymmetric theories at non-zero temperature, a topic of interest in 
recent developments [p8fl. 
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Figure Captions 

Figure 1: Mean lattice results for the free-energy density of hot SU(3) 
gauge theory from Ref . [[J . Here Si att refers to the free-energy divided by the 
free-energy of an ideal gas of gluons. 

Figure 2: The divergent perturbative free-energy density of SU (3) gauge 
theory given in Eq.@. Starting from the lowest curve at (a/27r) ' 5 = 0.24, 
one has iV = 2, 5, 3, 4. 

Figure (3a): The resummed perturbative free-energy density of hot SU (3) 
gauge theory for N = 3, 4 and 5, using a one-loop running coupling, the ref- 
erence values Cq = 1, Xq = 1, and the renormalization scale c = 1. The 
curves move upwards as iV increases. 

Figure (3b): The fifth order resummed perturbative free-energy density 
of Fig. (3a) at three different renormalization scales, c = 0.5, 1 and 2. The 
free-energy density increases with increasing c. 

Figure 4: The one-loop (upper curve) and approximate two-loop running 
couplings for SU(3) gauge theory at the renormalization scale c = 1. 

Figure (5a): The resummed perturbative free-energy density of hot SU (3) 
gauge theory for N = 3, 4 and 5, using a two-loop running coupling, the ref- 
erence values Co = 1, xq = 3, and the renormalization scale c = 1. The 
curves move upwards as N increases from 3 to 5. 

Figure (5b): The fifth order resummed perturbative free-energy density 
of Fig. (5a) at three different renormalization scales, c = 0.5, 1 and 2. The 
free-energy density increases with increasing c. 

Figure 6: The [1,2], [2,1] and [2,2] Borel-Pade' approximants to the 
perturbative free-energy density, with a two-loop running coupling, and the 
renormalization scale c = 1. Starting with the lowest curve at x = 5, one has 
[2,1], [1,2], [2,2]. 

Figure 7: The fifth order resummed perturbative free-energy denisty of 
SU(3) gauge theory coupled to Nf fermions as discussed in Sec. (7). Starting 
from below at x — 5, the curves label Nf = 6, 5, 4, 3, 0, 2, 1. 



24 



Figure (8a): The divergent perturbative free-energy density of QED, given 
in Sec. (9). Starting from the lowest curve and moving upwards, one has 
N = 2, 5, 3, 4. 

Figure (8b): The resummed perturbative free-energy density of QED. 
The curves move upwards as N increases from 3 to 5. 
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Fig. 6 
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